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THE  DEVELOPMENT  OF  RAPID  PREDICTIVE  METHODS  FOR 
THREE-DIMENSIONAL  TRANSONIC  FLOW  FIELDS 
ABOUT  FIGHTER  BOMBER  AIRCRAFT 


Anthony  J.  Crisalli,  Stephen  S.  Stahara, 
Jack  N.  Nielsen,  and  John  R.  Spreiter 


1.  INTRODUCTION 


The  safe  and  controllable  separation  of  external  stores,  either 
powered  or  unpowered,  from  modern  fighter-bomber  aircraft  remains  a 
current  problem  of  significant  military  importance  for  the  U.S.  Air 
Force.  Operational  limitations  due  to  dangerous  store  release  character¬ 
istics  can  reduce  aircraft  mission  effectiveness  and  survivability. 
Consequently,  development  of  accurate  engineering  predictive  methods  for 
describing  store  trajectories  fills  a  critical  dual  role.  Use  of  such 
methods  not  only  enhances  the  performance  and  safety  of  weapon  delivery, 
but  also  provides  a  means  for  reducing  the  time  required  for  both  wind- 
tunnel  and  full-scale  flight  tests  relating  to  store  certification 


prog  rams , 


Recognition  of  this  need  has  lead  Nielsen  Engineering  6  Research, 
Inc.  (NEAR)  through  AFFDL  support  to  develop  the  most  comprehensive 
engineering  predictive  method  currently  available  for  determining  store 
separation  characteristics  at  subcritical  subsonic  speeds.  A  similar 
program  for  developing  the  analogous  capability  at  supersonic  speeds  is 
currently  underway.  At  transonic  speeds,  however,  the  problem  becomes 
significantly  more  complicated  because  of  the  essential  nonlinear 
character  of  the  flow.  Nevertheless,  the  need  for  an  accurate  store 
trajectory  predictive  technique  at  transonic  speeds  is  even  more  com¬ 
pelling  as  unusual  store  separation  phenomena  are  more  likely  to  occur  in 
this  speed  regime. 


This  report  presents  an  overall  account  of  the  progress  made  on  an 
essential  component  of  the  transonic  store  trajectory  problem:  the 
prediction  of  transonic  flow  fields  about  wing-body/pylon  combinations 
representative  of  modern  fighter-bomber  aircraft.  The  emphasis  here  is 
on  the  development  of  predictive  techniques  which  are  sufficiently 
economical  from  the  standpoint  of  computer  cost  to  enable  an  aircraft 


designer  to  employ  the  techniques  over  a  wide  range  of  Mach  number, 
angle  of  attack,  and  geometric  parameters. 

The  report  presents  several  candidate  transonic  predictive  methods 
and  evaluates  them  by  comparing  their  predictions  of  various  flow  field 
quantities  with  experimental  data.  The  experimental  data  are  taken  from 
a  parallel  wind-tunnel  test  phase  of  the  project,  designed  specifically 
for  evaluating  the  predictive  methods.  The  experimental  data  have  been 
documented  and  archiv?d  in  several  test  reports  {ref«.  l^’O .  The  exten¬ 
sive  comparisons  of  the  theoretical  predictions  from  the  various 
predictive  methods  with  the  data  comprise  the  major  portion  of  this 
report. 

In  the  following  sections,  work  on  the  transonic  flow  field  pre¬ 
diction  problem  as  performed  under  APOSR  Contract  No.  F4<1620-75-C-0047 
will  be  reviewed  in  its  entirety,  in  such  a  way  as  to  highlight  the 
evolution  of  the  various  predictive  methods  from  the  initiation  of  the 
research  project  to  the  present.  Both  experimental  and  theoretical 
phases  are  covered.  Many  comparisons  of  the  predictive  methods  with 
experimental  data  are  shown  and  used  to  draw  conclusions  about  the  rela¬ 
tive  merits  of  the  various  methods.  A  subsequent  report  will  deal  with 
the  application  of  the  flow  field  predictive  techniques  discussed  in 
this  report  to  the  computation  of  external  store  loadings  at  transonic 
speeds. 


2.  OVERVIEW 

In  this  section  we  present  an  overview  of  the  theoretical  and  experi¬ 
mental  phases  of  the  project  insofar  as  they  relate  to  the  computation  of 
three-dimensional  transonic  flow  fields  about  transonic  wing-body  con¬ 
figurations  with  pylons  under  the  wing  or  fuselage. 

The  goal  of  the  theoretical  phase  is  to  develop  a  transonic  flow 
field  predictive  technique  which  is  sufficiently  economical  and  accurate 
to  be  used  as  a  basis  for  computing  store  loadings  and  trajectories.  In 
particular,  the  force  and  moment  distributions  on  a  store  which  is 
separated  from  the  parent  aircraft  (l.e.  the  wing-body/pylon  combina¬ 
tion)  can  be  found  by  first  calculating  the  flow  field  in  the  vicinity 
of  the  parent  aircraft  without  the  store  in  the  flow  field.  The  force 
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and  moment  distribution  on  the  store  can  then  be  determined  by  inputing 
the  velocity  components  at  the  position  the  store  occupies  in  the  flow 
into  a  store  loading  calculation.  Thus,  the  ability  to  predict  the  flow 
field  in  the  vicinity  of  the  wlng-body/pylon  combination,  without  store, 
is  the  essential  first  step  in  the  computation  of  store  loadings  and 
trajectories. 

In  view  of  the  above,  the  theoretical  phase  of  the  project  has 
focused  on  the  development  and  testing  of  several  flow-field  predictive 
methods  about  wing-body/pylon  configurations  and,  more  basically,  wing- 
body  configurations.  These  flow-field  predictive  methods  are  now  listed 
chronologically  in  the  order  in  which  they  have  been  investigated  in  the 
project: 

1.  The  classical  transonic  equivalence  rule 

2.  An  extended  transonic  equivalence  rule 

3.  Three-dimensional  linear  solutions  based  on  paneling  methods 

4.  A  three-dimensional  transonic  smal l-dlsturbance  finite-dif¬ 
ference  method 

In  this  report  application  of  the  above  methods  to  the  flow-field 
predictive  problem  are  described,  compared  against  experimental  data, 
and  evaluated. 

The  purpose  of  the  parallel  experimental  phase  of  the  project  is  to 
provide  data  for  critically  analyzing  and  evaluating  the  various  flow- 
field  predictive  methods.  In  two  4T  wind-tunnel  experiments  (ref.  1  and 
refs.  2-4),  an  idealized  P-16  model  was  tested  over  a  range  of  angles  of 
attack  at  Mach  numbers  in  the  transonic  regime.  Flow-field  data  as  well 
as  surface  pressure  and  force/moment  distributions  were  obtained.  This 
report  deals  exclusively  with  the  flow-field  data  from  these 
experiments . 

In  the  following  section  the  experimental  and  theoretical  phases  of 
the  project  will  be  presented  in  detail.  In  the  next  section,  section  3, 
a  description  of  the  wind-tunnel  tests  and  data  obtained  is  presented. 

In  sections  4  through  9,  analyses  and  evaluations  of  the  predictive  tech¬ 
niques  discussed  above  are  given.  Section  10  contains  an  analysis  of  the 
effects  of  wind-tunnel  wall  interference  on  the  experimental  data. 
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rinally,  section  11  presents  the  overall  conclusions  and  possible  future 
directions  suggested  by  the  Investigation  thus  far. 


3.  THE  EXPERIMENTAL  DATA 

The  experimental  data  to  be  used  in  the  evaluation  of  the  transonic 
flow-field  predictive  methods  is  taken  from  two  AEOC  4T  wind  tunnel 
tests.  These  tests  are  described  in  detail  in  the  data  reports  (refs.  1 
and  2-4).  Here  we  give  a  brief  description  of  the  experiments  with 
emphasis  on  the  flow-field  data.  The  question  of  the  validity  of  the 
experimental  data  for  the  evaluation  of  theoretical  methods  is  also 
briefly  discussed. 

3.1  Description  of  Experiments 

Both  experiments  utilized  the  same  idealized  F-16  wing-body  model 
in  configurations  with  and  without  pylons.  (Some  configurations 
contained  stores,  but  that  data,  which  involves  store  loadings,  will  be 
utilized  in  the  subsequent  report.)  Figure  1  shows  the  body-fixed 
coordinate  system  which  is  used  throughout  this  report.  Figure  2  shows 
planform  and  cross-sectional  views,  including  key  geometrical  locations 
of  the  wing-body/pylon/store  model.  Figure  3  shows  a  photograph  of  the 
wing-body  model  in  the  4T  tunnel  and  the  conical  flow-field  survey  probe 
which  was  used  to  obtain  flow-field  data  in  the  vicinity  of  the 
wing-body  model. 

In  the  first  transonic  wind-tunnel  experiment  (ref.  1),  the  wing- 
body  model  was  employed  in  the  AEDC  4T  tunnel  without  pylons.  The  pri¬ 
mary  purpose  of  this  experiment  was  to  obtain  transonic  flow-field  data, 
including  axial  velocity,  upwash  velocity,  and  sidewash  velocity, 
throughout  the  region  surrounding  the  wing-body  model.  The  data  were 
obtained  at  very  closely-spaced  intervals  in  order  to  capture  the  steep 
variations  in  flow-field  quantities  which  occur  in  the  vicinity  of  shock 
waves.  Both  inner  and  outer  flow-field  surveys  were  obtained.  The 
inner  flow-field  data  were  obtained  to  evaluate  the  various  theoretical 
predictive  methods  discussed  in  this  report.  The  outer  flow-field  data, 
taken  on  a  cylindrical  surface  as  far  from  the  model  centerline  as 
possible,  were  obtained  to  provide  measured  outer  boundary  conditions  as 
input  to  the  theoretical  model  in  order  to  evaluate  wind-tunnel  inter- 


ference  effects.  The  locations  of  the  Inner  and  outer  flow-field  grids 
are  shown  in  figures  4  and  5,  respectively.  In  this  report,  a  repre¬ 
sentative  subset  of  the  flow-field  surveys  from  the  first  test  are 
employed.  The  selected  data  involves  two  angles  of  attack,  a  >  0°,  5^, 
and  the  three  transonic  free-stream  Mach  numbers,  M^  ■  0.925,  0.975, 
1.025. 

A  thorough  survey  of  these  4T  experimental  results  of  reference  1 
has  verified  that  the  test  parameters  were  exceptionally  well  selected 
in  view  of  providing  as  wide  a  range  of  transonic  phenomena  as  possible. 
The  data  display  flow  conditions  from  subcrltical  to  slightly  super¬ 
critical,  to  strongly  supercritical,  to  mildly  supersonic,  as  were 
desired.  For  the  two  subsonic  free-stream  Mach  numbers,  figure  6 
illustrates  this  fact  and  displays  the  growth  of  the  supersonic  pockets 
on  the  pressure  and  suction  sides  of  the  wing.  The  results  are  for  an 
(x,z)  plane  located  at  the  spanwise  location  y  >  2  inches  (25%  semispan) 
and  are  for  the  wing-body  model  having  a  thickness-only  wing  profile 
based  on  the  actual  thickness  distribution  of  the  F-16  wing.  The  figure 
on  the  top  indicates  the  extent  of  the  supersonic  zone  at  M  •  0.925  for 
the  three  angles  of  attack,  while  corresponding  results  for  M^  ■  0.975 
are  shown  in  the  bottom  plot.  The  symbol  M^^  denotes  the  local  Mach 
number.  Since  the  vertical  limits  of  the  inner  flow  surveys  was 
1  <  Izl  <  5  inches,  the  maximum  lateral  locations  of  the  larger  super¬ 
sonic  pockets  on  the  suction  side  were  beyond  the  last  inner  survey 
location  at  z  ■  5  inches.  However,  only  for  the  M^  ■  0.975,  a  •  5° 
case  did  the  pocket  extend  out  to  the  outer  flow  survey  location  at  z  ■ 
14  Inches.  These  results  indicate  the  extreme  sensitivity  of  the  flow 
at  supercritical  conditions.  Analogous  results  for  M^  ■  1.025  are  shown 
in  figure  7  which  displays  the  variation  and  growth  of  the  embedded 
subsonic  pocket.  Of  particular  note  in  both  figures  6  and  7  is  that, 
at  modest  angles  of  attack,  flow  conditions  on  the  pressure  side  of  the 
wing  remain  primarily  subsonic  for  a  wide  range  of  conditions. 

In  the  second  experiment  (refs.  2-4),  the  wing-body  in  combination 
with  various  pylon  and  store  arrangements  was  tested  at  several  Mach 
numbers  and  angles  of  attack  in  the  AEDC  4T  tunnel.  The  purpose  of  this 
experiment  was  to  obtain  systematic  inner  and  outer  flow-field  data  as 
well  as  store  surface  pressures,  forces  and  moments.  The  wing-body/ 
pylon/store  configurations  used  in  the  flow-field  testing  are  shown  in 


figure  8.  The  symbol  x  In  figure  8  denotes  the  position  in  the  cross 
plane  of  the  flow-field  survey.  The  flow-field  data  includes  axial 
velocity,  upwash,  and  sldewash  velocities.  This  report  utilizes  only 
flow-field  data  in  the  vicinity  of  the  following  configurations  (con¬ 
figurations  21,  22,  and  26  in  fig.  8): 

1.  Wing-body 

2.  Wing-body  with  a  single  pylon  on  each  wing 

3.  Wing-body  with  pylon  on  fuselage 

This  report  presents,  for  the  purpose  of  theoretical  evaluation,  a 
representative  subset  of  the  flow-field  data  for  the  above  configura¬ 
tions  from  the  second  experiment  (refs.  2-4).  The  selected  data 
Includes  two  angles  of  attack,  a  ■  0°,  5°,  and  transonic  free-stream 
Mach  numbers,  M_  ■  0.95,  1.05. 

3.2  Validity  of  the  Experimental  Data 


The  use  of  experimental  data  to  evaluate  predictive  theories  must 
be  qualified  by  several  considerations  pertaining  to  the  quality  of  the 
data.  These  are: 

1.  Data  uncertainties  due  to  instrumentation  errors 

2.  Tunnel  flow-field  quality  imperfections 

3.  Flow-field  lateral  asymmetry 

4.  Wind-tunnel  wall  interference  effects 

In  references  1  and  2,  it  was  seen  that  the  probable  maximum  uncer¬ 
tainties  due  to  the  combined  effect  of  the  first  three  considerations 
above  are  ♦  .015  for  axial  velocity  as  normalized  by  the  free-stream 
velocity,  *  .40°  for  upwash  angle,  and  *  .40°  for  sldewash  angle. 
Therefore,  the  maximum  precision  that  one  can  reasonably  demand  from  a 
theoretical  prediction  method  is  that  it  agree  with  the  experimental  data 
to  within  these  uncertainty  tolerances. 


The  importance  of  assessing  the  fourth  consideration,  wind-tunnel 
wall  interference  effects,  stems  from  the  fact  that  wall  interference 
effects  have  been  assumed  to  be  negligible  in  the  theoretical  predic¬ 
tions  whereas  for  some  of  the  data,  particularly  that  taken  at  Mach 
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numbers  close  to  one,  they  may  not  be.  The  detailed  discussion  of  inter* 
ference  effects,  which  relates  to  some  of  the  computational  results  of 
the  following  sections  as  well  as  the  data  from  a  companion  wind-tunnel 
test  in  the  AEDC  16T  tunnel  (ref.  1)  is  deferred  to  section  10.  At  this 
point  we  simply  limit  ourselves  to  the  conclusions  regarding  wall  inter¬ 
ference  effects  which  are  drawn  in  section  10: 

1.  For  <  0.975,  the  experimental  data  are  relatively  free  of 
wall  interference  effects. 

2.  For  ^  1.025,  wall  interference  effects  are  not  significant 
over  most  of  the  wing.  However,  they  may  be  significant  before 
the  wing  leading  edge  and  in  the  vicinity  of  the  wing  trailing 
edge. 


4.  THE  CLASSICAL  TRANSONIC  EQUIVALENCE  RULE 

4.1  Description  of  Classical  Transonic  Equivalence  Rule 

The  classical  transonic  equivalence  rule  can  be  expressed 
mathematically  as 


2  #  Q  2  #  u 


“  *2,B 


where  «  is  the  composite  perturbation  potential  for  the  transonic  flow 
field,  and  ^2  t'  ^2  B  solutions  of  the  two-dimensional 

Laplace  equation  as  shotm  in  figure  9.  Here,  is  the  two-dlmen- 

sional  incompressible  flow  solution  for  the  translation  of  the  cross 
section  (lifting  problem),  ^  corresponding  solution  for  expansion 
(or  contraction)  of  the  cross  section  (thickness  problem),  and  ^2  b 
solution  for  the  expansion  (or  contraction)  of  a  circular  cylinder  having 
identical  cross-sectional  area.  Finally,  is  the  solution  to  the  non¬ 
linear  transonic  small-disturbance  equation  for  axisymmetric  flow  past 
the  equivalent  body  of  revolution. 
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4.2  Comparisons  of  Classical  Transonic  Equivalence  Rule  with  Data 


Extensive  comparisons  of  the  classical  transonic  equivalence  rule 
with  the  experimental  data  obtained  in  the  initial  4T  test  (ref.  1)  have 
been  made.  Only  a  sample  of  those  results  will  be  presented  in  this 
report.  Representative  examples  of  these  comparisons  are  shown  in 
figures  10  and  11.  Table  1  contains  an  index  of  these  comparisons  and 
the  corresponding  page  numbers. 


Figure  10  shows  the  axial  velocity,  upwash  and  sidewash  comparisons 
at  ■  0.925,  a  ■  0°  at  the  longitudinal  survey  whose  (y.z)  coordinates 
are  (y.z)  -  (2.-1).  Figure  11  shows  the  same  comparisons  for  ■  0.975 
Cross-sectional  planform  sketches  of  the  wing-body  showing  the  approxi¬ 
mate  position  of  this  survey  are  shown  below  (all  dimensions  are  in 
inches) : 


Radius  *1.33 


urvey 


On  the  basis  of  the  comparisons  indicated,  as  well  as  numerous  other 
comparisons  not  shown,  the  following  observations  can  be  drawn: 


Axial  velocity:  The  classical  equivalence  rule  prediction 
is  satisfactory  for  points  ahead  of  the  leading  edge  of  the  wing 
tip.  i.e..  X  ■  16.1  inches.  Downstream  of  the  leading  edge  of 
the  wing  tip  the  predicted  axial  velocity  shows  large  disagree¬ 
ment  with  the  data. 


Upwash:  The  classical  equivalence  rule  predictions  are 

generally  satisfactory  except  for  a  tendency  in  some  cases  to 
overpredict  the  experimental  data  just  before  the  wing 
trailing  edge,  as  illustrated  in  figure  11(b). 


10.0  / 
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Sidewash:  The  classical  equivalence  rule  predictions  are  satis 

factory  at  locations  ahead  of  the  leading  edge  of  the  wing  tip 
(X  <■  16.1  inches).  Beyond  this  point,  the  prediction  is 
generally  unsatisfactory,  exhibiting  a  severe  overprediction 
and/or  underprediction  of  the  data. 


4.3  Shortcomings  of  Classical  Transonic  Equivalence  Rule 


The  overall  conclusions  that  can  be  drawn  from  the  above  compari¬ 
sons,  as  well  as  many  others  (not  shown)  made  over  a  broad  range  of 
Mach  numbers,  angles  of  attack,  and  survey  locations  is  that  for  the 
configuration  geometry  considered,  the  classical  transonic  equivalence 
rule  is  reliable  at  locations  upstream  of  the  leading  edge  of  the  wing 
tip  (x  •>  16.1  inches),  but  generally  unreliable  downstream  of  the 
leading  edge  of  the  wing  tip  in  its  predictions  of  axial  velocity,  up- 
wash,  and  sidewash. 


Initially,  it  was  believed  that  the  poor  performance  of  the  equiva¬ 
lence  rule  beyond  x  ■  16.1  inches  could  be  ascribed  to  numerical 
problems  caused  by:  (1)  the  sharp  break  in  the  wing  planform  at  the 
wing  tip,  (2)  trai 1 ing-edge  singularities  caused  by  inaccuracies  in  the 
wing  profile  second  derivatives,  and  (3)  the  presence  of  a  boundary 
layer  of  signficant  thickness  in  the  vicinity  of  the  wing  trailing  edge 


In  an  attempt  to  correct  these  problems,  the  geometry  of  the  wing 
was  smoothed  as  shown  in  the  illustration  below. 


Wing-tip  V 

smoothing:  ^ 

(planform  view) 


Trailing  edge  V.«, 
smoothing : 

(wing  profile  view) 


where  the  solid  line  represents  the  actual  wing  and  the  dotted  line 
indicates  the  smoothing  used.  Additionally,  the  wing  profile  second 
derivatives  were  smoothed  using  a  least-squares  cubic  spline  represents 
tion  of  the  wing  profile. 


The  above  modifications  ameliorated  the  problem  only  slightly.  Fur¬ 
ther  analysis  Indicated  that  the  disagreement  between  theory  and  data 
resulted  from  a  more  fundamental  cause.  Consider  the  unsatisfactory 
agreement  of  the  axial  velocity  u/V^  as  predicted  by  the  original  tran¬ 
sonic  equivalence  rule  with  data  for  a  nonlifting  flow  at  ■  0.975  as 

shown  in  figure  10(a).  For  a  -  0«  u/V^  can  be  written  as 

JL  ^  ^ 

V  V  V 

The  variation  of  Ug/V^  is  shown  in  figure  12(a)  while  that  for  ~ 

u^  is  shown  in  figure  12(b).  It  is  apparent  that  Ug/V^  is  quite 

smooth  and  behaves  qualitatively  like  the  data.  However,  ~ 

U2  does  not  and  exhibits  a  discontinuous  derivative  at  x  *  16.1 

inches  (the  location  of  the  wing-tip  leading  edge).  An  examination  of 
the  analytic  expression  for  U2  ^./V  -  U2  g/V^  shows  that  this  is  caused 

by  discontinuities  in  the  wing  planform  (i.e.,  at  the  wing  tip  and  the 
trailing  edge)  which  induce  true  three-dimensional  effects — the  predic¬ 
tion  of  which  are  beyond  the  capability  of  slender-body  theory  which  forms 
the  basis  of  the  classical  equivalence  rule.  The  reason  behind  this 
breakdown  is  that  slender-body  theory  performs  satisfactorily  as  long  as 
the  wing-span  is  increasing  and  the  planform  is  slender.  These  restric¬ 
tions  become  even  more  important  when  a  ^  0,  since  slender-body  theory 
predicts  that  beyond  the  point  of  maximum  span,  a  trailing  vortex  sheet 
emanates  from  the  wing  and  extends  downstream  parallel  to  the  flow  direc¬ 
tion.  Beyond  that  point  of  maximum  span,  the  axial  lift  distribution 
is  zero.  Consequently,  the  problem  is  not  numerical  but  theoretical, 
and  can  be  resolved  only  by  modifying  or  extending  the  classical 
equivalence  rule  so  as  to  overcome  this  shortcoming  of  slender-body 
theory.  A  procedure  for  accomplishing  this  will  be  discussed  in  the 
next  section. 


5.  THE  EXTENDED  TRANSONIC  EQUIVALENCE  RULE 

5.1  Description  of  Extended  Transonic  Equivalence  Rule 

The  extended  transonic  equivalence  rule  can  be  represented  by  the 
following  formula: 
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♦  “  *  *8  '  *3,8 

Here  a,  ^  is  the  solution  to  the  three-dimensional  linear  flow 
^^3,  t+a 

around  the  configuration  at  an  arbitrary  angle  of  attack  a,  4'^  is  the 
usual  solution  to  the  ax i symme tr ic  nonlinear  transonic  equation  around 
the  equivalent  body,  at  zero  angle  of  attack,  and  g  is  the  axisym- 
metric  linear  solution  to  the  flow  around  the  equivalent  body  at  zero 
angle  of  attack. 

The  rationale  behind  the  extended  equivalence  rule  will  now  be  dis¬ 
cussed.  In  section  4,  it  was  seen  that  the  classical  equivalence  rule 
does  not  properly  account  for  planform  discontinuities  at  the  wing  tip 
and  trailing  edge.  This  is  due  to  the  fact  that  the  components  ^ 
and  g  of  the  classical  equivalence  rule  are  two-dimensional  solutions 
computed  in  crossflow  planes  perpendicular  to  the  longitudinal  axis  of 
the  aircraft.  In  order  to  represent  the  proper  three-dimensional  be¬ 
havior  of  the  solution  in  these  regions,  it  therefore  seems  appropriate 
to  replace  ^2  8  their  three-dimensional  counterparts 

♦j^^^and  g,  respectively.  This  yields  the  extended  equivalence 
rule.  It  should  be  noted  that  the  composite  velocity  potential  ^  in 
the  above  equation  is  asymptotically  equivalent  to  ^g  at  Infinity  and 
asymptotically  equivalent  to  on  the  aircraft  surface. 

An  advantage  of  the  extended  equivalence  rule  is  the  fact  that  the 
solution  method  for  the  component  4-,  is  already  available  in  the 

form  of  operational  codes  (refs.  5  and  6)  based  on  panel  methods.*  A 
conceivable  disadvantage  of  the  extended  equivalence  rule  is  the  fact 
that  the  1  inear- theory  solution,  t*a'  calculable  at  free- 

stream  Mach  numbers  near  1.  Extensive  calculations  which  are  discussed 
below,  however,  have  shown  that  the  linear  theory  paneling  methods  can 
be  used  to  provide  useful  results  even  at  free-stream  Mach  numbers  quite 
close  to  1. 


*  The  supersonic  paneling  code  of  reference  6  has  been  modified  slightly 
for  this  report.  The  code  of  reference  6  includes  special  treatment 
of  nonlinearities  due  to  shocks.  That  treatment  has  been  removed  for 
the  work  presented  in  this  report. 
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In  the  following  subsections  the  comparisons  of  the  extended  tran¬ 
sonic  equivalence  rule  with  data  from  reference  1  are  presented.  The 
notation  used  is  that  the  extended  transonic  equivalence  rule  predictions 
are  represented  by  the  solid  symbol  with  tick  mark,  the  experimental 

data  by  the  open  symbol  O,  and  the  1  inear-theory  predictions  by  the 
solid  symbol  #.  The  1  inear- theory  predictions  will  be  discussed 
separately  in  section  6. 

5.2  Comparison  of  Extended  Transonic  Equivalence  Rule  with 
Wing-Body  Data,  ■  0.925 

Typical  comparisons  between  the  extended  equivalence  rule  and  the 
wing-body  alone  data  at  ■  0.925,  a  ■  0°  are  shown  in  figures  13  and 
14.  In  these  figures  the  predicted  axial  velocity  upwash,  and  sidewash 
are  compared  with  data  at  surveys  located  at  (y,z)  ■  (2,-1)  and  (y,z)  > 
(4,-1).  The  approximate  positions  of  these  surveys  are  shown  in  the 
sketch  below. 


Radius  > 1. 33 


It  is  seen  in  figures  13  and  14  that  the  agreement  between  the  extended 
equivalence  rule  and  the  data  is  almost  exact  for  these  cases.  Many 
other  comparisons  (not  shown  in  this  report)  at  a  •  0°  and  •  5®  and  over 
a  range  of  various  survey  locations  confirm  that  the  agreement  between 
the  extended  equivalence  rule  flow-field  predictions  and  experimental 
data  for  this  configuration  is  generally  very  good  at  M„  -  0.925.  Thus, 
at  this  Mach  number,  the  extended  equivalence  rule  provides  a  distinct 
improvement  over  the  classical  transonic  equivalence  rule. 

5.3  Comparison  of  Extended  Transonic  Equivalence  Rule  with 
Wing-Body  Data,  M^  ■  0.975 

We  now  consider  the  comparisons  of  the  extended  equivalence  rule 
with  experimental  data  at  M„  >  0.975.  Comparisons  for  angles  of  attack 
a  ■  0®,  5®  at  survey  locations  (y,z)  ■  (0,-2),  (2,-1),  (2,-3),  (4,-1) 
are  provided  in  figures  15  to  22. 
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A  Sketch  of  the  wing-body  illustrating  the  approximate  (y,2)  loca¬ 
tions  noted  above  Is  provided  below. 


On  the  basis  of  the  comparisons  indicated  in  these  figures,  the 
following  general  conclusions  regarding  the  extended  equivalence  rule 
flow-field  predictions  at  ■  0.975  can  be  made. 

1.  Axial  velocity:  The  extended  equivalence  rule  predictions  are 
satisfactory  at  axial  locations  up  to  roughly  the  leading  edge 
of  the  wing  tip.  Beyond  this  location  predictions  are  quali¬ 
tatively  correct  but  generally  overpredict  the  data. 

2.  Upwash:  The  extended  equivalence  rule  predictions  are  satis¬ 
factory  upstream  of  the  trailing  edge  of  the  wing  except  just 
under  the  fuselage  {y,z)  ■  (0,-2)  where  the  comparison  is  poor 
downstream  of  the  leading  edge  of  the  wing  tip. 

3.  Sidewash:  The  extended  equivalence  rule  predictions  are  satis¬ 
factory  upstream  of  the  leading  edge  of  the  wing  tip.  Beyond 
this  point  the  theory  overpredicts  and  then  underpredicts  the 
da  ta . 

Comparison  of  the  extended  equivalence  flow  field  predictions  with 
the  classical  equivalence  rule  flow-field  predictions  for  the  cases 
shotm  provides  the  conclusion  that  the  extended  rule  is  definitely 
superior. 

5.4  Comparison  of  Extended  Transonic  Equivalence  Rule  with 
Wing-Body  Data,  ■  1.025 

Analogous  comparisons  of  the  extended  equivalence  rule  and  data  are 
presented  in  figures  23  to  30  for  angles  of  attack  a  *  0^,  5°  at  the 
surface  locations  (y,2)  ■  (0,-2),  (2,-1),  (2,-3),  (4,-1).  The  following 
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general  conclusions  regarding  the  extended  equivalence  rule  flow-field 
predictions  at  ■  1.025  can  be  drawn  on  the  basis  of  the  above  com¬ 
parisons: 

1.  Axial  velocity:  Extended  equivalence  rule  predictions  are  satis¬ 
factory  at  axial  locations  upstream  of  the  leading  edge  of  the 
wing  tip.  Beyond  this  location  they  greatly  overpredict  the 
data.  This  however,  may  be  due  to  wall  interference  effects 
which  are  discussed  in  section  10. 

2.  Upwash:  Extended  equivalence  rule  predictions  are  satisfactory 
upstream  of  the  wing  trailing  edge  except  for  the  survey  just 
under  the  fuselage  (y,z)  -  (0,-2)  where  beyond  the  wing  tip 
leading  edge  the  theory  underpredicts  and  then  overpredicts 

the  data. 

3.  Sidewash:  Extended  equivalence  rule  predictions  are  generally 
satisfactory  up  to  the  wing  trailing  edge  except  for  the  survey 
at  about  midspan  (y,z)  ■  (4,-1).  For  this  case  at  a  ■  5^  [see 
fig.  30(c) 1,  the  theory  greatly  overpredicts  the  data  at  the 
wing  leading  edge. 

Overall,  the  extended  equivalence  rule  is  generally  superior  to  the 
classical  equivalence  rule  for  the  mildly  supersonic  case  ■  1.025. 

5.5  Analysis  of  Extended  Equiva  *.ice  Rule  Behavior 
As  stated  previously,  the  extended  transonic  equivalence  rule  is 

*  "  *3,t+a  ^  *B  ■  *3,B 

Since  is  the  full  three-dimensional  linear  solution  for  the  wii.g- 

body,  it  follows  that  ♦g  -  ♦3^3  represents  a  nonlinear  correction  to 
that  linear  solution.  The  magnitude  of  the  nonlinear  correction  term  is 
represented  in  figures  13  to  30  by  the  differences  between  the  symbols 
^  and  9. 

With  regard  to  the  extended  equivalence  rule  axial  velocity  pre¬ 
dictions  at  ■  0.975  we  observe  that  between  the  leading  edge  of  the 
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wing  tip  and  the  wing  trailing  edge  the  nonlinear  correction  term 

♦n  ~  ♦in  is  too  large  near  the  body  fsee  figs.  15(a),  16(a),  19(a), 20(a), 

D  J  I  B 

and  21(a)].  It  should  be  noted  that  the  cross-sectional  area  of  the 
equivalent  body  changes  40%  in  this  region.  It  appears  that  such  a 
sudden  change  in  area  is  too  much  for  the  equivalence  rule  to  handle  at 
subsonic  speeds  so  close  to  ■  1. 

At  ■  1.025,  the  tendency  of  the  extended  equivalence  rule  to  over¬ 
predict  the  data  with  regard  to  axial  velocity  is  due  to  an  entirely 
different  cause.  From  figure  23  to  30,  we  note  that  the  nonlinear  cor¬ 
rection  term  is  fairly  small  so  that  the  disagreement  with  the  data  is 
primarily  due  to  the  linear  component  rather  than  the  nonlinear  correc¬ 
tion.  The  problem  with  the  linear  theory  which  occurs  at  low  supersonic 
free-stream  Mach  numbers  is  related  to  multiple  reflections  and  oscilla¬ 
tions  occur  ing  in  the  panel  method  (refs.  5  and  6)  and  can  be  corrected 
by  modifications  of  the  existing  panel  technique. 

Regarding  the  extended  equivalence  rule  upwash  predictions  at 
subsonic  or  supersonic  speeds,  we  note  that  the  predictions  are  quite 
reasonable  for  most  of  the  surveys  shown.  Under  the  fuselage  the  com¬ 
parison  between  theory  and  data  for  upwash  is  poor  between  the  leading 
edge  of  the  wing  tip  and  the  wing  trailing  edge.  Again,  the  reason  for 
this  is  that  the  nonlinear  correction  term  is  too  large.  Regarding  the 
extended  equivalence  rule  sidewash  predictions,  we  note  that  under  the 
wing  (e.g.,  at  (y,z)  ••  (2,-1)  and  (4,-1)  the  sidewash  predictions  are 
poor  between  the  leading  edge  of  the  wing  tip  and  the  wing  trailing  edge, 
again  because  of  overcorrection  by  the  nonlinear  transonic  term. 

5.6  Overall  Conclusions  Regarding  Extended  Transonic 

Equivalence  Rule 

The  overall  conclusions  drawn  from  these  numerous  comparisons  is  that 
the  extended  equivalence  rule  is  generally  superior  to  the  classical 
equivalence  rule  for  the  prediction  of  flow-field  quantities.  Compari¬ 
sons  of  the  theoretical  predictions  with  data  exhibit  almost  exact 
agreement  at  -  0.925.  However,  discrepancies  with  the  extended  equiv¬ 
alence  rule  theory  exist  at  0.975  and  >  1.025  and  seem  to  be  due 

primarily  to  overcorrection  by  the  nonlinear  transonic  term  ♦g  -  g. 
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at  *  0.925.  Additional  comparisons  (not  shown  in  this  report)  at 
a  a  0°,  5°  and  over  a  range  of  survey  locations  confirm  the  following 
conclusions ; 

1.  The  linear  theory  and  the  extended  equivalence  rule  are 
virtually  identical  in  their  predictions  of  axial  velocity, 
upwash  and  sidewash  at  -  0.925. 

2.  The  comparisons  of  the  linear-theory  flow-field  Quantities  with 
experimental  data  show  good  agreement  at  ■  0.925. 

6.3  Comparison  of  Linear  Theory  with  Wing-Body  Data, 

-  0.950 

Typical  comparisons  between  the  linear  theory  and  the  experimental 
data  at  ■  0.95,  a  ■  0°,  5°  are  shown  in  figures  31  to  39.  In  these 
figures  the  predicted  axial  velocity,  upwash  and  sidewash  are  compared 
with  data  at  the  survey  locations  (y,z)  ■  (0,-2.208),  (0,-2.958), 
(.750,-2.208),  (3.5,-1.230),  (3.5,-1.980),  and  (4.25,-1.230).  The 
approximate  (y,z)  positions  are  designated  by  the  symbol  "x"  in  the 
sketch  below. 


In  the  above  sketch  the  dotted  lines  indicate  the  positions  which  the 
pylons  and  stores  would  occupy  if  they  were  present.  (In  section  7,  the 
data  with  pylons  present  will  be  given,  and  then  the  above  data  will  be 
useful  in  assessing  the  effect  of  the  pylon  on  the  flow  field).  The 
following  conclusions  are  drawn  from  these  comparisons: 

1.  Axial  velocity:  The  comparisons  of  linear  theory  with  data 
are  satisfactory.  Some  discrepancies  can  be  seen  for  most 
cases  between  the  leading  edge  of  the  wing  tip  (x  ■  16.1  inches) 
and  the  wing  trailing  edge. 
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2.  Upwash:  The  conparisons  of  linear  theory  with  data  are  satis¬ 
factory.  The  theory  does  overpredict  the  data  near  the  wing 
leading  and  trailing  edges  in  some  cases. 

3.  Sidewash:  The  comparisons  of  linear  theory  with  data  are 
satisfoctory.  The  theory  does  overpredict  the  data  near  the 
wing  leading  edge  and  the  leading  edge  of  the  wing  tip. 

6.4  Comparison  of  Linear  Theory  with  Wing-Body  Data, 

-  0.975 

We  now  discuss  the  conparisons  of  the  linear  theory  with  experi¬ 
mental  data  for  the  wing-body  at  >  0.975.  Figures  15  to  22  show 
comparisons  of  axial  velocity,  upwash  and  sidewash  for  angles  of  attack 
o  ■  0°  and  5°  at  the  survey  locations  (y,2)  -  (0,-2),  (2,-1),  (2,-3), 
(4,-1).  A  sketch  of  the  approximate  location  of  these  positions  was 
shown  in  section  5.3.  The  following  conclusions  can  be  drawn  from  these 
comparisons: 

1.  Axial  velocity:  Linear  theory  predictions  are  satisfactory  up 
to  roughly  the  leading  edge  of  the  wing  tip.  Beyond  this  point 
the  linear  prediction  is  unsatisfactory  and  generally  under¬ 
predicts  the  data. 

2.  Upwash:  Linear  theory  predictions  are  satisfactory  except 
downstream  of  the  wing  trailing  edge.  In  some  cases  the  theory 
significantly  overpredicts  the  data  at  the  wing  leading  edge.* 

3.  Sidewash:  Lineary  theory  predictions  are  satisfactory  in  all 
cases  except  for  some  disagreement  near  the  wing  leading  edge 
for  the  (y,z)  ■  (4,-1)  case. 

A  comparison  of  the  above  conclusions  with  the  corresponding 
conclusions  for  the  extended  equivalence  rule  for  ■  0.975  presented  in 


*The  discrepancy  between  the  theory  and  data  for  the  upwash  at  the 
leading  edge  of  the  wing  appears  to  be  due  to  an  insufficient  number 
of  chordwise  lifting  panels  on  the  wing. 
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section  5.3  shows  that  the  linear  theory  is  superior  to  the  extended 
equivalence  rule  in  the  prediction  of  upwash  and  sidewash  and  somewhat 
inferior  in  the  prediction  of  axial  velocity. 


6.5  Comparison  of  Linear  Theory  with  Wing-Body  Data, 

-  1.025 

The  comparisons  of  the  linear  theory  with  experimental  data  for  the 
wing-body  alone  at  >  1.025  are  included  in  figures  23  to  30.  It  is 
seen  in  these  figures  that  the  predictions  sometimes  display  small  oscil 
latory  behavior  at  axial  locations  ahead  of  the  leading  edge  of  the 
wing  tip  (X  -  16.1  inches)  and  large  oscillations  at  locations  beyond 
that  point.  It  is  believed  this  behavior  is  caused,  in  part,  by  multipl 
wave  reflections  occurring  in  the  supersonic  paneling  method  used  in 
reference  6  for  M  near  1. 

The  following  conclusions  can  be  drawn  from  these  comparisons: 

1.  Axial  velocity:  Linear- theory  predictions  are  satisfactory  up 
to  the  leading  edge  of  the  wing  tip.  Beyond  this  point  they 
greatly  overpredict  the  data. 

2.  Upwash:  The  1  inear-theory  predictions  are  satisfactory  up  to 
the  wing  trailing  edge. 

3.  Sidewash:  The  1  inear- theory  predictions  are  satisfactory  up  to 
the  wing  trailing  edge  except  for  some  disagreement  near  the 
wing  leading  edge  for  the  (y,z)  >  (4,-1)  case. 

Comparing  these  conclusions  with  those  for  the  extended  transonic 
equivalence  rule  in  section  5.4,  it  is  seen  that  under  the  wing  there  is 
little  difference  between  the  upwash  and  sidewash  predictions  of  the  two 
theories  at  ■  1.025.  However,  under  the  fuselage  the  linear  theory 
yields  superior  upwash  predictions. 
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6.6  Comparison  of  Linear  Theory  with  Wing-Body  Data, 

-  1.050 

Typical  comparisons  between  linear  theory  and  the  experimental 
data  at  ■  1.05,  a  -  0°,  5°  are  shown  in  figures  40  to  48.  In  these 
figures  the  predicted  axial  velocity,  upwash,  and  sidewash  are  compared 
with  data  at  the  same  survey  locations  as  in  section  6.3  (see  sketch  in 
section  6.3).  The  following  conclusions  are  drawn  from  these  comparisons 

1.  Axial  velocity:  Li  near- theory  tends  to  overpredict  the 
deceleration  of  the  flow  under  the  model.  The  predictions 
also  tend  to  give  a  much  stronger  compression  wave  downstream 
of  the  wing  trailing  edge  than  is  seen  in  the  data. 

2.  Upwash:  The  1  inear- theory  predictions  are  satisfactory  up  to 
the  wing  trailing  edge. 

3.  Sidewash:  The  1  inear- theory  predictions  are  satisfactory  up  to 
the  wing  trailing  edge.  The  theory  does  tend  to  overpredict 
the  data  near  the  wing  leading  and  trailing  edges. 

It  is  seen  from  figures  40  to  48  that  some  minor  oscillatory  behavior 
is  still  present  in  the  predictions  although  it  is  less  than  for  the  case 
-  1.025  discussed  in  section  6.5.  However,  some  of  the  oscillations 
are  not  small  (see  figs.  41(a),  42(b),  43(c),  44(c),  45(c),  and  48(c)]. 

6.7  Overall  Conclusions  Regarding  Linear  Theory 

On  the  basis  of  the  above  comparisons  it  is  seen  that  the  linear 
theory  yields  satisfactory  upwash  and  sidewash  predictions  over  a  range 
of  survey  positions  under  the  fuselage  and  wing  at  free-stream  Mach 
numbers  relatively  close  to  one.  The  axial  velocity  as  predicted  by  the 
linear  theory  is  reasonably  good  at  M^  ■  0.925  and  M^  ■  1.05,  but 
increasingly  deteriorates  beyond  the  wing  tip  leading  edge  as  M^  -*  1. 
Furthermore,  for  M^  >  1,  the  linear  theory  predictions  of  all  three 
velocity  components  show  an  increasingly  oscillatory  behavior  as  M^  ■*  1 
from  the  supersonic  side. 


22 


When  compared  with  the  performance  of  the  classical  transonic  equiv¬ 
alence  rule  and  the  extended  transonic  equivalence  rule,  it  is  seen  that 
the  linear  theory  yields  more  reliable  predictions  of  both  upwash  and 
sidewash  over  a  broad  range  of  conditions  for  0®  <  a  <  5®  and  <  0.975 
M  >  1.025.  It  is  also  seen  that  the  axial  velocity  predictions  of  the 
linear  theory  are  better  than  those  of  the  classical  equivalence  rule 
but  not  so  good  as  those  of  the  extended  rule.  However,  it  appears  that 
axial  velocities  may  have  only  a  second-order  effect  on  the  normal  and 
sidewash  forces  and  the  pitching  and  yawing  moments  of  a  separated  store 
Thus,  on  balance,  the  linear  theory  is  considered  the  best  of  the  flow- 
field  predictive  methods  in  the  low  angle  of  attack  range,  <  0.975  or 
M  >  1.025,  for  the  wing-body  configuration  considered  for  angles  of 
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attack  up  to  5®. 

It  should  be  noted  that  the  above  conclusions  are  based  on  compari¬ 
sons  with  experimental  data  for  wing-body  configurations  without  pylons. 
In  the  next  section  the  linear  theory  will  be  compared  with  wing-body/ 
pylon  experimental  data. 

7.  COMPARISON  OF  THREE-DIMENSIONAL  LINEAR  THEORY  PREDICTIONS 
WITH  WING-BODY/PYLON  DATA 

7.1  Description  of  the  Wlng-Body/Pylon  Configuration 
and  Survey  Locations 

In  this  section  we  examine  the  effects  on  the  flow  field  of  adding 
wing  and  fuselage  pylons  to  the  wing-body.  As  discussed  in  section  3, 
the  second  4T  experiment  Involved  the  following  configurations:  wing- 
body,  wing-body  with  t%«i  wing  pylons,  and  wing-body  with  fuselage  pylon. 
The  survey  (y,z)  positions  associated  with  the  wing-body/pylon  configure 
tions  are  described  below.  The  comparisons,  with  corresponding  page 
numbers,  are  contained  in  Table  3. 

For  the  wing-body  with  two  wing  pylons  configuration,  the  survey 
{y,z)  positions  used  in  the  comparisons  are  (y,z)  ■  (3.5,-1.230),  (3.5, 
-1.980),  (4.250,-1.230).  The  approximate  (y,z)  positions  are  designated 
by  the  symbol  "X"  in  the  sketch  below. 
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In  the  above  and  following  sketch,  the  dotted  lines  represent  the 
positions  which  the  stores  would  occupy  if  they  were  present.  Note  that 
the  surveys  are  located  along  the  centerlines  of  stores  for  which 
loading  data  were  obtained  in  the  experiment. 

Note  that  the  above  survey  positions  are  among  the  survey  positions 
for  wing-body  data  discussed  in  sections  6.3  and  6.6.  These  wing-body 
data  together  with  wl ng-body/pylon  data  allow  the  effect  on  the  flow 
field  of  adding  the  pylon  to  be  assessed  experimentally. 

For  the  wing-body  with  fuselage  pylon  configuration,  the  survey 
(y,z)  positions  used  in  the  graphical  comparisons  are  (y,z)  ■  (0-2.208), 
(0.-2.958),  (.750,-2.208).  The  approximate  (y,z)  positions  are  shown  in 
the  sketch  below. 


Note  again  that  these  survey  positions  are  also  among  the  survey 
positions  for  the  wing-body  configurations  discussed  in  sections  6.3  and 


The  1 Inear- theory  predictions  are  compared  with  experimental  data 
involving  the  above  wl ng-body/pylon  configurations  and  designated 
survey  locations  in  figures  49  to  63,  as  indexed  in  Table  3.  Two  free- 
stream  Mach  numbers,  M^  •  0.95,  1.05,  and  two  angles  of  attack,  o  *  0°, 
5®,  are  treated.  The  effect  of  adding  the  pylon  is  treated  in  section  8. 
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7.2  Discussion  of  Linear-Theory  Predictions  for  Wlng- 
Body/Pylon  Configurations 

A  scan  of  ail  the  comparisons  discussed  above  supports  the  follow¬ 
ing  conclusions: 

1.  At  M  ■  0,950  (figs.  49-54  and  61-63),  the  upwash  and  sldewash, 

<l>  * 

as  computed  by  the  linear  theory,  are  In  satisfactory  agreement 
with  the  data.  The  axial  velocity  shows  good  agreement  with 
the  data  upstream  of  the  wing  tip-wing  pylon  trailing  edge 
region.  However,  the  linear  theory  does  not  model  correctly 
the  flow  expansion  and  contraction  in  this  region  due  to  the 
presence  of  an  embedded  pocket  of  supersonic  flow. 

2.  For  -  1.050  (figs.  55-60),  the  1 Inear-theory  predictions  show 
varying  degrees  of  oscillatory  behavior  in  axial  velocity,  up¬ 
wash,  and  sldewash  when  the  pylons  are  present.  The  oscillations 
generally  span  the  data.  The  oscillating  behavior  is  expecially 
pronounced  beyond  x  ■  16.5  inches.  Comparing  these  wing-body/ 
pylon  predictions  with  the  corresponding  predictions  for  the 
wing-body  at  the  same  Mach  number  as  discussed  in  section  6.6, 

it  is  seen  that  the  presence  of  the  pylon  has  intensified  the 
oscillatory  behavior.  As  mentioned  before,  these  oscillations 
are  associated  with  the  use  of  panel  methods  near  M^  ■  1. 

7.3  Overall  Conclusions  Regarding  Linear  Theory 

On  the  basis  of  the  results  shown  in  this  section,  it  is  concluded 
that  for  wing-body/pylon  configurations  at  M^  ■  0.95,  the  linear  theory 
predicts  upwash  and  sldewash  reasonabi y  wel 1 .  Axial  velocity  is  not 
satisfactorily  predicted  beyond  x  ■  16.5  Inches.  At  M^  •  1.050  the 
1  inear- theory  predictions  are  generally  unsatisfactory  because  of  oscil¬ 
lations  in  the  axial  velocity,  upwash,  and  sldewash  predictions.  Without 
the  pylons  present,  the  oscillations  are  much  milder.  Possible  improve¬ 
ments  to  the  linear  theory  in  this  respect  will  be  discussed  in 
section  11. 


8.  EFFECT  OF  PYLON  ADDITION  ON  THE  TRANSONIC  FLOW  FIELD 

Figures  64  to  78,  as  indexed  in  Table  4,  show  the  direct  effect  of 
pylon  addition  to  the  flow  field.  In  each  of  these  figures,  the  symbol 
O  represents  the  differences  (WBP(EXP)  -  WB(EXP)1  and  the  symbol  X 
represents  the  difference  tWBP(LT)  -  WB(LT)),  where 

WBP(EXP)  -  experimentally  measured  flow-field  quantity  (as 
specified  on  the  y-axis)  in  the  vicinity  of  the 
wl ng-body/pylon  combination. 

WB(EXP)  -  exper imen tal 1 y  measured  flow-field  quantity  in  the 
vicinity  of  the  wing-body. 

WBP(LT)  -  linear  theory  prediction  of  flow  field  quantity  in  the 
the  vicinity  of  the  wing-body/pylon  combination. 

WB(LT)  -  linear  theory  prediction  of  flow  field  quantity  in  the 
vicinity  of  the  wing-body. 

Thus,  the  difference  (WBP(EXP)  -  WB(EXP)J  directly  represents  the 
experimentally-measured  effect  of  pylon  addition  on  a  particular  flow 
field  quantity.  The  difference  (WBP(LT)  -  WB(LT)1  similarly  represents 
the  effect  of  pylon  addition  as  predicted  by  3-D  linear  theory. 

These  comparisons  provide  numerous  physical  insights  into  the  effect 
of  pylon  addition  on  the  flow  field.  For  example,  examining  figure  64(b) 
which  illustrates  the  effect  of  pylon  addition  on  upwash,  we  observe 
that  the  leading  edge  of  the  pylon  turns  the  flow  sharply  downward  while 
the  trailing  edge  of  the  pylon  turns  the  flow  sharply  upward. 

A  thorough  examination  of  these  graphs  provides  the  following  general 
conclusions: 

1.  At  M  •  0.950  (figs.  64-69  and  76-78),  the  effect  of  pylon  addi¬ 
tion  on  upwash  and  sidewash  is  satisfactorily  predicted  by  3-D 
linear  theory.  The  effect  of  pylon  addition  on  axial  velocity 
is  not  always  satisfactorily  predicted  by  the  linear  theory  in 
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tlhre  code  is  limited  to  subsonic  free-stream  Mach  numbers.  It  can  treat 
only  clean  wing-body  configurations,  i.e.,  without  pylons  or  stores. 

9.2  Utility  of  Three-Dimensional  Nonlinear  Transonic 

Theory  Resul ts 

There  are  several  reasons  for  comparing  the  prediction  of  the  three- 
dimensional  transonic  code  with  experiment. 

1.  The  predictions  are  interesting  in  their  own  right.  The  compari¬ 
sons  to  be  discussed  below  are  probably  the  first  extensive  comparisons 
of  three-dimensional  transonic  small-disturbance  results  with  experi¬ 
mental  flow-field  data. 

2.  Since  the  three-dimensional  transonic  code  can  be  expected  to  be 
more  accurate  than  all  the  rapid  predictive  methods,  its  prediction  when 
compared  with  data  fixes  an  error  bracket  which  rapid  predictive  methods 
cannot  be  expected  to  surpass. 

3.  The  three-dimensional  transonic  code  provides  a  tool  to  analyze 
wind-tunnel  wall  interference  effects.  This  aspect  will  be  discussed  in 
section  10. 

4.  The  three-dimensional  transonic  code  can  serve  as  a  basis  upon 
which  to  develop  improved  rapid  predictive  methods.  This  will  be 
discussed  more  fully  in  section  11. 

9.3  Comparison  of  Three-Dimensional  Nonlinear  Transonic 
Theory  Result  with  Data 

The  three-dimensional  transonic  code  was  utilized  to  obtain  flow 
field  predictions  in  the  vicinity  of  the  wing-body  configuration  for  the 
following  cases: 

1.  -  0.950,  a  -  0°,  5® 

2.  M_  -  0.975,  o  -  0°,  5® 

In  each  of  the  above  cases,  axial  velocity,  upwash,  and  sidewash 
predictions  were  obtained  at  a  variety  of  survey  locations  in  the 
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vicinity  of  the  wing-body,  including  all  survey  locations  contained  in 
the  data  reports  (refs.  1-4). 

In  this  report,  the  three-dimensional  transonic  code  predictions  will 
be  compared  with  experimental  data  only  for  the  cases  ■  0.975,  a  ■  0°, 
5°.  Although  predictions  at  many  survey  locations  were  obtained,  only 
the  representative  surveys  (y,z)  ■  (0,-2),  (2,-1),  (2,-3),  (4,-1),  (4,-3), 
(7,-1),  and  (7,-3)  will  be  used  for  the  axial  velocity,  upwash  and  side- 
wash  comparisons.  These  comparisons  are  shown  in  figures  79  to  92.  A 
detailed  listing  of  the  comparisons  shown  in  this  report  with  corres¬ 
ponding  page  numbers  is  contained  in  Table  5.  The  approximate  location 
of  these  surveys  is  shown  in  the  sketch  below. 
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Based  on  these  comparisons  the  following  conclusions  can  be  drawn 
regarding  capabilities  of  the  three-dimensional  transonic  small-dlstrlbu- 
tion  result: 

1.  The  predictions  of  axial  velocity,  upwash  and  sidewash  compare 
satisfactorily  in  most  cases  with  the  data  except  at  the  wing  trailing 
edge. 


2.  The  only  predictions  which  are  quantl ta tlvel y  only  fair  (but, 
nevertheless,  qualitatively  correct)  are  the  upwash  prediction  at  (y,z)  • 
(0,-2)  and  a  •  5®  (fig.  06(b) 1  and  the  sidewash  prediction  at  (y,z)  ■ 
(2,-1)  (figs.  80(c)  and  87(c)).  The  probable  cause  of  these 
discrepancies  is  that  these  surveys,  which  are  located  close  to  the 
fuselage,  are  more  strongly  affected  by  the  approximation  used  in  the 
present  version  of  the  code  to  treat  the  fuselage  boundary  conditions, 
i.e,  that  of  transferring  the  boundary  conditions  from  the  actual 
circular  fuselage  to  a  square  box  fuselage.  In  fact,  linear  theory  is 
somewhat  better  than  transonic  theory  at  these  points  as  can  be  seen  by 
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conparing  figures  16(c)  and  80(c),  19(b)  and  86(b),  and  20(c)  and  87(c). 
Comparisons  (not  shown)  of  the  3-D  transonic  theory  with  data  for  surveys 
farther  away  from  the  fuselage,  (y,2)  ■  (0,-3),  (0,-4),  (0,-5),  show 
excellent  agreement. 

3.  The  predictions  of  upwash  and  sidewash  by  the  three-dimensional 
transonic  code  are  comparable  in  accuracy  to  the  predictions  of  upwash 
and  sidewash  by  the  linear  theory.  This  supports  the  belief  that 

for  the  angle  of  attack  and  wing-body  geometry  considered  here,  non¬ 
linear  transonic  effects  are  small  insofar  as  upwash  and  sidewash  are 
concerned.  Therefore,  it  is  appropriate  to  use  a  (rapid)  linear  theory 
procedure  to  predict  these  quantities  at  low  angles  of  attack  up  to 
-  0.975. 

4.  The  prediction  of  axial  velocity  by  the  three-dimensional  tran¬ 
sonic  code  is  a  great  Improvement  In  accuracy  over  the  linear  theory 
prediction  beyond  the  point  x  ■  16.5  inches.  This  is  the  region  of  the 
shock  where  nonlinear  transonic  effects  can  be  expected  to  be  strongest. 
It  appears  inappropriate  to  use  a  1 Inear- theory  code  to  predict  axial 
velocity  in  the  shock  region  of  a  transonic  flow  field,  as  would  be 

ex pec  ted . 

10.  ASSESSMENT  OP  WIND-TUNNEL,  WALL-INTERFERENCE  EFFECTS 

The  assessment  of  wind-tunnel,  wall-interference  effects  In  the 
transonic  experimental  data  is  important  because  all  predictive  methods 
discussed  in  this  report  use  a  free-air  outer  boundary  condition.  Con¬ 
sequently,  if  the  effect  of  the  wind  tunnel  walls  on  the  data  is  not 
negligible,  there  could  be  a  signf leant  discrepancy  between  the  theo¬ 
retical  predictions  and  the  data. 

The  fact  that  the  three-dimensional  nonlinear  transonic  results 
(which  use  a  free-air  boundary  condition)  was  observed  in  section  9  to 
agree  well  with  the  experimental  data,  supports  the  conclusion  that  the 
experimental  data  at  M_  £  0.975  are  relatively  free  of  wall  interference 
effects.  What  differences  do  exist  between  theory  and  experiment  may  be 
due  to  minor  wall-interference  effects,  wind-tunnel  flow  irregularities, 
measurement  inaccuracies,  model  fabrication,  viscous  effects,  or  approxi¬ 
mations  in  the  theory. 
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The  assessment  of  wall-interference  effects  at  M  ■  1.025  must  pro- 
ceed  differently  than  for  ■  0.975  because  the  current  version  of  the 
Bailey-Ballhaus  code  cannot  handle  supersonic  free-stream  Mach  numbers. 
Instead  the  method  developed  under  this  contract  and  reported  in  ref¬ 
erence  9  is  used.  This  method  employs  measured  outer  boundary  conditions, 
as  obtained  in  the  wind-tunnel  tests  of  references  1-4,  as  the  outer 
boundary  condition  for  the  transonic  axlsymmetric  equivalent-body  calcula¬ 
tion.  The  measured  outer  boundary  condition  results  for  the  equivalent 
body  are  then  compared  with  corresponding  results  using  a  free-air  outer 
boundary  condition.  This  comparison  is  shown  for  equivalent  body  surface 
pressure  in  figure  93  for  M^  ■  1.025.  It  is  seen  in  figure  93  that  the 
only  regions  of  disagreement  are  in  the  vicinities  of  x  ■  11.0  inches 
(i.e.,  before  the  wing  leading  edge  over  most  of  the  span)  and  x  >  17.0 
inches.  Therefore  this  lends  support  to  the  conclusion  that  wall  inter¬ 
ference  effects  in  the  case  M  <■  1.025  are  not  significant  over  the  bulk 
of  the  wing  but  may  be  significant  only  in  front  of  the  wing  leading  edge 
and  behind  x  -  17.0  inches.  This  result  may  apply  more  to  u/V^  than  to 
downwash  and  sidewash  angles. 

Finally,  we  note  that  above  conclusions  regarding  the  relatively 
minor  Influence  of  wind-tunnel  wall  interference  at  M  -  0.975  and  1.025 
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are  consistent  with  the  conclusions  drawn  in  reference  10.  There,  the 
conclusion  that  wall-interference  effects  had  only  a  minor  influence  on 
the  flow  field  was  drawn  on  the  basis  of  a  comparison  between  4T  experi¬ 
mental  data  and  16T  experimental  data. 

11.  OVERALL  CONCLUSIONS  AND  FUTURE  DIRECTIONS 

The  main  conclusion  to  emerge  from  the  investigation  described  in 
this  report  is  that  the  three-dimensional  linear  theory  predictive 
technique  yields  generally  satisfactory  rapid  predictions  of  upwash  and 
sidewash  in  a  transonic  flow  field  around  the  wing-body  or  wing-body/ 
pylon  conf igurations  considered,  over  a  range  of  transonic  free-stream 
Mach  numbers  and  angles  of  attack. 

For  the  subsonic  free-stream  Mach  numbers  considered  (M^  ■  0.925, 
0.950,  and  0.975)  and  the  angles  of  attack  considered  (o  ■  0®  and  5°),  the 
linear  theory  yielded  satisfactory  agreement  in  upwash  and  sidewash  up¬ 
stream  of  the  wing  trailing  edge  for  both  the  wing-body  and  wing-body/ 


fields.  The  three-dimensional  transonic  code,  which  does  incorporate 
nonlinear  effects,  was  seen  to  yield  good  agreement  with  data  in  the 
vicinity  of  the  wing-body  alone  combination  at  ■  0.975,  a  ■  0°,  5°. 

As  presently  constituted,  it  cannot  handle  configurations  involving 
pylons.  Additionally,  it  is  not  computationally  economical  and  is 
presently  limited  to  subsonic  free-stream  Mach  numbers.  Nevertheless, 
the  possibility  exits  that  the  three-dimensional  transonic  code  could  be 
incorporated  on  a  limited  basis  as  part  of  a  larger  scheme  for  the 
computation  of  transonic  flow  fields.  Such  a  scheme  can  be  formulated 
as  follows: 

WBP  -  (WBP)lt  * 

Here 

WBP  ■  predicted  flow  field  quantity  (axial  velocity,  upwash, 
or  sidewash)  in  the  vicinity  of  the  wing-body/pylon 
conf iguration. 

(WBP)j^^  ■  flow  field  quantity  as  predicted  by  the  three-dimensional 
linear  theory  for  the  wing-body/pylon  configuration 

(WB)nl  ■  flow  field  quantity  as  predicted  by  the  three-dimensional 
transonic  theory  for  the  wing-body  configuration 

(WB)lt  *  flow  field  quantity  as  predicted  by  the  three-dimensional 
linear  theory  for  the  wing-body  configuration. 

The  above  scheme  handles  wing-body/pylon  combinations  and  it 
utilizes  the  three-dimensional  transonic  code  (no  pylons)  to  produce  non¬ 
linear  transonic  effects.  Thus  the  extension  of  the  3-D  transonic  method 
to  handle  pylons  would  be  circumvented  by  this  device.  The  problem  of 
economy  in  the  above  formulation  can  be  solved  by  utilizing  perturbation 
techniques  (refs.  11-12)  which  enable  one  to  approximate  an  entire  range 
of  3-D  transonic  predictions  from  just  two  or  three  3-D  transonic  pre¬ 
dictions.  For  example,  the  3-D  transonic  code  could  be  used  to  predict 
the  flow  field  for  the  cases  a  >  0°,  5°  and  then,  from  these  two  cases, 
an  entire  range  of  nonlinear  predictions  could  be  obtained  via  a  highly 
economical  and  accurate  perturbation  technique.  The  accuracy  of  this 
procedure  could  be  checked  by  comparing  the  procedure's  flow  field  pre- 
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dictions  at  a  ■  10°  or  15°,  or  by  obtaining  wind-tunnel  test  data  at 
a  m  10°  and  15°  and  using  the  data  as  a  check.  Finally,  the  extension 
of  the  existing  3-D  transonic  code  to  supersonic  free-stream  Mach  numbers 
requires  both  the  implementation  of  a  supersonic  free-air  outer  boundary 
conditions  as  well  as  some  grid  modifications  and  should  be  carried  out. 
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TABLE  2.-  COMPARISONS  OF  EXTENDED  TRANSONIC  EOL’IVALENCE  RULE 
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TABLE  3.-  COMPARISONS  OF  LINEAR  THEORY  WITH  DATA  FOR 
REFERENCES  2-4 


TABLE  4.-  COMPARISONS  OF  EFFECT  OF  PYLON  ADDITON  AS  GIVEN  BY 
LINEAR  THEORY  WITH  EFFECT  OF  PYLON  ADDITION 
AS  GIVEN  BY  DATA  OF  REFERENCES  2-4 


TABLE  5.-  COMPARISONS  OF  3-D  TRANSONIC  PREDICTIONS 
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Figure  4,-  General  grid  for  inner  flow-field  measurements 
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re  10.-  Comparison  of  measured  and  theoretical  flow  field  quantities 
r  wing  of  wing-body  combination  at  0.925,  a  «  0®,  y-2  in.,  z  - -1  ii 


(b)  Upwash  angle,  deg. 


Figure  13.-  Comparison  of  measured  and  theoretical  flow  field  quantities 
for  wing-body  combination:  a  ■  0°,  ■  0.925.  y  •  2  in.,  z  ■  -1  in. 


Model  station,  in. 


Figure  15.-  Comparison  of  measured  and  theoretical  flow  field  quantities 
for  wing-body  combination;  a  ■  o'’,  M  "  0.975,  y  ■  0  in.,  z  »  -2  in. 
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Figure  19.-  Concluded 
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gure  32.-  Continued 


Figure  33.-  Comparison  of  measured  flow  field  quantities  with  prediction 
of  thre'i-dimensional  linear  theory  for  wing-body  combination, 
a  -  0°.  -  0.95,  y  -  0.750  in.,  z  -  -2.208  in. 


(c)  Sidewash  angle,  deg. 
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Figure  35.-  Concluded. 


(b)  Up%>rash  angle,  deg 


Figure  39.-  Comparison  of  measured  flow  field  quantities  with  prediction 
of  three-dimensional  linear  theory  for  wing-body  combination, 
a  -  5°,  M  -  0.95,  y  -  4.250  in.,  z  -  -1.230  in. 
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Figure  43,-  Comparison  of  measured  flow  field  quantities  with  prediction 
of  three-dimensional  linear  theory  for  wing-body  combination, 
a  ■  0^.  M  ■  1,-5.  y  “  3,5  in,,  z  ■  -1,230  in. 


Figure  49.-  Comparison  of  measured  flow  field  quantities  with  prediction 
of  three-dimensional  linear  theory  for  wing-body  combination  with  wing 
pylons,  a  -  0°,  M  ■  0.95,  y  •  3.5  in.,  z  ■  -1.230  in. 
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re  51.-  Comparison  of  measured  flow  field  quantities  with  prediction 
three-dimensional  linear  theory  for  wing-body  combination  with  wing 
pylons,  a  »  0°,  -  0.95,  y  »  4.250  in.,  z  -  -1.230  in. 


re  52.-  Comparison  of  measured  flow  field  quantities  with  prediction 
three-dimensional  linear  theory  for  wing-body  combination  with  wing 
pylons,  a  -  5°,  M,  ■  0.95,  y  -  3. 5  in.,  z  »  -1.230  in. 


Figure  52.-  Continued 
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Figure  54.-  Continued 
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(a)  Axial  velocity  ratio. 

Figure  59.-  Comparison  of  measured  flow  field  quantities  with  prediction 
of  three-dimensional  linear  theory  for  wing-body  combination  with  wing 
pylons,  a  ■  5  ®,  *  1.05,  y  ■  3.5  in.,  z  ■  -1.980  in. 
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Figure  68.-  Continued 
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Figure  79.-  Concluded. 


(b)  Upwash  angle,  deg. 


Figure  03.-  Concluded 


Piqur«  86.-  Coroparison  of  measured  flow  field  quantities  for  wing-body 
combination  with  predictions  of  Bailey-Ballhaus  transonic  code, 
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Figur®  06.-  Concluded 
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(b)  Upwash  angle,  deg 
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(b)  Upwash  angle,  deg 


SVMBOLS 


FP 

LE 

LT 

NL 


PLE 

PTE 

TE 

u,v,w 

u/V 

m 

v_ 

WB 

WB/FP 

WP 

WB/WP 


x,y,z 

xt,yt,2t 


pressure  coefficient  on  surface  of  axisymmetric 
equivalent  body 

fuselage  pylon 

wing  leading  edge 

linear  theory  predictive  n»ethod  (i.e.,  3-D  panel 
method,  refs.  5-6) 

free-stream  Mach  number 

nonlinear  theory  predictive  method  (i.e.,  3-D  transonic 
code,  refs.  7-8) 

pylon  leading  edge 

pylon  trailing  edge 

wing  trailing  edge 

local  velocity  components,  positive  along  the  positive 
x,y,z  directions,  respectively 

axial  velocity  ratio 

free-stream  velocity 

wing-body  configuration 

configuration  consisting  of  wing-body  with  fuselage 
pylon 

wing  pylon 

configuration  consisting  of  wing-body  with  pylon  on  each 
wing 

body  fixed  Cartesian  coordinate  system  shown  in  figure  1 

tunnel-fixed  Cartesian  coordinate  system  directed 
parallel  and  perpendicular  to  the  tunnel  centerline. 

The  origin  is  coincident  with  the  aircraft  model  nose 
at  zero  angle  of  attac)c.  Positive  directions  are 
identical  to  the  body  axis  system. 

model  angle  of  attac)i,  deg.,  as  shown  in  figure  1 

upwash  angle,  tan~^  (w/u) ,  deg. 

sidewash  angle,  tan~^  (v/u) ,  deg. 

difference  between  wing-body /pylon  and  wing-body  axial 
velocities 
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difference  between  winq-body/pylon  and  wing-body 
upwashes 

difference  between  wing-body/pylon  and  wing-body  side- 
washes 

small  perturbation  velocity  potential  for  various 
theories  considered  in  report 

small-perturbation-velocity-potential  components  of 
classical  transonic  equivalence  rule  described  in 
section  4 

small-perturbation-velocity-potential  components  of 
extended  transonic  equivalence  rule  described  in 
section  S 
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